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CHAPTER  1 


I N  T  H  OIMJC  T  1 O  N 


Thig  report  presents  an  analysis  of  a  series  of  radiation  measurements  of  an 
underground  ammunition  magazine  converted  into  a  shelter.  Data  obtained  by  the 
U.S.  Naval  Radiological  Defense  Laboratory  were  compared  with  previously  pub¬ 
lished  iheoretical  and  experimental  work, and  the  protection  afforded  by  the  shelter 
if  it  were  exposed  to  a  full  infinite  licld  of  contamination  of  either  fallout  or  cobalt- 
60  was  estimated.  (  ) 


BACKGROUND 

The  subject  of  this  analysis,  the  Naval  Radiological  Defense  Laboratory  (NRDL) 
skelter  and  the  experimemal  measurements  performed  upon  it,  is  described  in  the 
paper,  "Preliminary  Measurements  of  Shielding  Effectiveness  of  an  Underground 
Fallout  Shelter."1  The  shelter  is,  in  effect,  a  quonset  hut  25  ft  wide  and  about 
49  ft  long,  buried  so  that  the  floor  is  nearly  12  ft  below  the  ground.  Earth  has  been 
packed  flat  over  the  roof  so  that  the  minimum  roof  cover  is  2  ft  10  in.  thick  at  the 
center  of  the  hut  arch. 

The  entrance  is  an  L-shaped  passageway  that  was  designed  to  have  no  line-of- 
sight  leakage  path  for  radiation.  It  is  of  circular  cross  section  7  ft  in  diameter. 

•ad  there  is  a  3/16-in.  thick  steel  door  covering  the  opening.  At  the  rear  of  the 
■kelter  is  a  vent  2  ft  in  diameter  with  its  center  3  ft  from  the  rear  wall,  which  is 
aaod  for  exhaust  and  for  an  emergency  escape  hatch.  Radiation  leakage  through 
tke  entrance  and  the  exhaust  vent  is  of  primary  interest  in  the  experiment. 

Measurements  of  the  radiation  doses  46  in.  above  the  floor  were  made  from 
Chre  simulated  source  areas:  (1)  an  annulus  extending  from  a  1  -  to  5-ft  radius 
Msrrotmding  the  shelter  vent,  (2)  an  annulus  extending  from  a  5  -  to  25- ft  radius 
also  surrounding  the  shelter  vent,  {3)  a  30  x  30-ft  area  over  the  roof,  (4)  a  25  x 
M-ft  area  centered  on  the  entranceway  with  the  door  covered,  and  (5).  the  same 


25  x  50-ft  area  centered  on  the  entranceway  but  with  the  door  cleared.  Area  sources 
were  simulated  by  circulating  a  130-curie  cobalt-60  source  at  a  uniform  rate  through 
2000  ft  of  polyethylene  tubing  laid  out  in  suitable  arrays.  The  radiation  penetrating 
into  the  shelter  was  measured  with  time- integrating  detectors. 


ORGANIZATION  OF  THIS  REPORT 

The  organization  of  this  report  is  as  follows.  Chapter  2  discusses  the  atten¬ 
uation  characteristics  of  fallout  and  cobalt  -60  radiation.  Chrptcr  3  presents  specific 
computations  for  extending  the  experimentally- measured  radiation  intensities  to 
infinite- field  conditions.  These  experimental  results  are  compared  with  previously 
published  theory.  Chapter  4  presents  the  conclusions  drawn  from  the  study. 
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CHAPTER  2 

COMPARISON  OF  FALLOUT  AND  COBALT-GO  RADIATION 

We  present  in  th.s  chapter  some  basic  information  on  the  attenuation  character¬ 
istics  of  fallout  and  cobalt -GO  radiation. 


FALLOUT  RADL4TION 

The  energy  distribution  of  fallout  contamination,  as  computed  by  Nelms  and 
2 

Cooper,  continuously  changes  with  time  after  burst.  The  present  convention  of 

using  the  1-hr  energy  distribution  in  computing  shelter  factors  was  proposed  by 
3 

Spencer,  who  reasoned  that  a  large  part  of  the  total  exposure  to  radiation  is  apt  to 
occur  during  the  first  few  hours  and  that  calculations  based  on  the  3pectrum  at  1  hr 
after  the  explosion  would  provide  a  conservative  estimate  of  the  penetrability. 

The  energy  distribution  of  fallout  at  1. 12  hr,  23.  8  hr,  and  9.  82  days  after 
explosion  is  illustrated  in  Figure  1.  Note  that,  while  the  energy  averages  about 
1. 7  Mev  at  early  and  late  limes,  a  connonent  at  about  0.  8  Mev  dominates  at  times 
on  the  order  of  one  day.  The  effect  of  this  energy  shift  is  illustrated  in  Figure  2 
(reproduced  from  ref.  3),  a  plot  of  the  relative  protection  provided  by  an  underground 
shelter  with  300  psf  of  earth  cover  vs  time  after  nuclear  explosion. 

The  penetration  of  radiation  into  a  structure  may  be  conveniently  described  by 
three  basic  geometric  arrangements,  as  shown  in  Figure  3.  First,  there  is  pene¬ 
tration  from  sources  adjacent  to  a  vertical  barrier  (Case  II);  second,  there  is  pene¬ 
tration  from  sources  lying  on  a  barrier  (Case  I);  and  third,  there  is  penetration 
from  sources  adjacent  to  a  horizontal  barrier  (Case  m).  The  first  two  geometric 
representations  may  be  treated  in  a  straightforward  manner,  since  most  of  the 
radiation  penetrating  the  barrier  comes  directly  from  the  sources  rather  than  from 
scattering  by  the  air.  The  third  geometric  representation,  the  penetration  of  radia¬ 
tion  from  sources  adjacent  to  a  horizontal  barrier,  also  depends  upon  the  fraction 
of  the  radiation  "available"  for  penetration  (direct  radiation  being  "unavailable"  for 

3 

penetration  in  this  case).  Spencer  has  found  that  an  infinite  plane  field  of  fallout 
oontaminatton  produces  approximately  9.3%  of  its  dose  as  skyshinc  at  3  ft  above  the 


RELATIVE  PROTECTION  FACTOR 


PHOTON  ENERGY  I  MEV) 


Figure  1.  Distribution  of  Dose  3  Ft  Above  Infinite  Plane  Contaminated  With 
Fallout  of  Various  Ages  (Volatile  Products  Removed) 
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Figure  2.  Variation  of  Protection  Factor  with  Time 
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Figure  3.  Penetration  of  1. 12-Hour  Fallout  and  Cobalt-60  Radiation 
Into  Concrete"  (from  Spencer3) 
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J  plane,  while  a  similar  field  of  coball-GO  contamination  produces  an  8.4%  skyshine 

j  component  at  3 -ft  height.  Skyshine  is  defined  here  as  air-scattered  radiation  in  the 

!  downward  direction. 
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COBALT -60  RADIATION 

i 

5 

The  cobalt -60  skyshine  fraction  has  been  verified  experimentally  by  Rexroad. 
With  a  cobalt -60  source  on  the  ground,  Rexroad  measured  the  dose  rates  at  heights 
of  1,  3,  and  6  ft  above  the  ground  for  various  source  distances.  His  analysis 
results  in  buildup  factors  of  1. 24,  1. 25,  and  1. 28  respectively  for  these  heights 
above  an  infinite  contaminated  plane.  Thus,  the  scattered  dose  rate  for  both  sky¬ 
shine  (radiation  entering  from  above)  and  ground-scattered  radiation  may  be  com¬ 
puted  to  be  19. 4%,  20%,  and  21. 9%  of  the  total  (direct  plus  scattered)  radiation. 

The  skyshine  portion  of  the  scattered  intensity  would  be  half  of  these  values  (sinee 
skyshine  radiation  may  only  penetrate  from  the  upper  2Jr  steradians  of  a  sphere 
surrounding  a  detector)  if  the  ground  were  assumed  to  be  made  up  of  air  compressed 
to  the  same  density  as  the  ground. 

Table  1  shows  there  is  good  agreement  between  Spencer's  theoretical  results 
and  Rexroad's  experimental  results.  Thus,  cobalt -60,  like  fallout  radiation,  creates 
a  skyshine  field  of  only  about  10%  of  the  total  dose  rate  at  3  ft  above  an  infinite  plane 
of  contamination.  Air-scattered  radiation  from  cobalt -60  may  therefore  be  used  as 
a  skyshine  simulant  if  its  attenuation  can  be  related  in  a  systematic  way  to  the  atten¬ 
uation  of  fallout  skyshine. 

TABLE  1 


SKYSHINE  FRACTION 


Height 

(ft) 

Computed  Fallout 
(Spencer*) 

Computed  Cobalt -60 
(Spencer*) 

Measured  Cobalt -60 
(Rexroadt) 

1 

9.8% 

8.6% 

9.7% 

3 

9.3% 

8.4% 

10.0% 

6 

8.9% 

8.3% 

10.9% 

♦See  Ref.  3. 
*SeeRef.  5. 
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The  attenuation  to  both  cobalt -GO  ami  fallout  radiation  afforded  by  a  concrete 
barrier  has  been  determined  by  Spencer  for  the  three  idealized  geometries,  as 
!  Figure  3  shows.  Because  the  ratio  of  the  attenuation  provided  by  a  barrier  to  both 
;  fallout  and  cobalt  radiation  varies  significantly  with  both  thickness  and  geometry, 

\  the  shelter  factors  measured  experimentally  with  cobalt  -60  may  be  corrected  for 
fallout  only  in  terms  of  these  three  individual  geometric  situations.  Since  the  varia- 
>  tion  is  regular  in  all  cases,  however,  the  use  of  an  appropriate  equivalent  thickness 
j  computation  will  allow  cobalt-measured  protection  factors  to  be  related  to  fallout 

r 

j  protection  factors. 

Table  2  summarizes  the  relative  attenuation  of  1. 12 -hr  fallout  and  cobalt-GO 
|  radiation  for  various  barrier  thicknesses.  The  ratios  show  that  results  obtained 
|  by  use  of  the  cobalt -60  simulant  will  generally  be  over -optimistic.  In  particular, 

[  since  most  of  the  radiation  reaching  the  underground  shelter  proceeds  by  the  Case 
j  m  mechanism  of  downward  scattering,  protection  factors  must  be  reduced  to  2/3 

j  of  the  observed  (cobalt)  results  to  yield  appropriate  values  for  1-hr  fallout  radiation. 


TABLE  2 


RELATIVE  ATTENUATION  OF  1.12-HOUR  FALLOUT  AND 
COBALT-60  RADIATION 


Barrier 

pef 

Case  I: 
on 

Contamination 
a  Barrier 

Case  II:  Contamination 
Next  to  Vertical  Barrier 

Case  III:  Contamination 
Next  to  Horizontal  Barrier 

FO 

Co-60 

Ratio 

FO  : 

Co-60 

Ratio 

FO 

Co-CO 

Ratio 

60 

.078 

.069 

1.13 

.260 

.275 

0. 94 

.067 

.046 

1.46 

100 

.020 

.0165 

1.21 

.082 

.082 

1.00 

.0082 

.0056 

1.46 

160 

.0057 

.0042 

1.34 

.0276 

.0235 

116 

.00106 

. 00071 

1.49 

200 

.00178 

. 00115 

1.50 

.0095 

.0069 

1.37 

. 00016 

— 

260 

300 

.00065 

.000185 

.00029 

.000076 

1.90 

2.45 

.0034 

.00122 

.00205 

. 00056 

1.66 

2.18 
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CHAPTER  3 


ANALYSIS 

This  chapter  presents  an  estimate  of  the  radiation  dose  that  would  have  been 
delivered  within  the  NRDL  shelter  if  it  had  been  exposed  to  infinite-field  contamina¬ 
tion.  Considerable  portions  of  this  analysis  are  based  directly  upon  the  experimen¬ 
tal  data  by  Shumway,  et  al.  ,*  obtained  by  exposing  the  shelter  to  a  simulated 
uniformly-contaminated  but  limited  field  of  cobalt  radiation.  No  attempt  has  been 
made  to  analyze  the  shelter  completely,  but  rather  attention  has  been  directed 
toward  the  critical  areas,  such  as  the  entranceway  and  the  vent. 

SHELTER  VENT  AREA 

The  roof  of  the  shelter  contains  a  vent  structure  protruding  from  the  shelter  to 
a  point  2.  5  ft  above  the  ground  surface.  This  vent  structure  is  composed  of  a  2-ft 
diameter  pipe  of  l/4-in.  wall  thickness  with  a  l/4-in.  iron  plate  cover.  Radiation 
penetration  of  the  vent  was  investigated  by  simulating  a  field  of  contamination  in  two 
annular  areas  centered  about  the  vent.  These  two  areas,  called  the  skyshlne  array 
and  the  lip  array,  were  of  25-ft  outer  radius  5-ft  inner  radius,  and  5-ft  outer  radius 
1-ft  inner  radius  respectively.  Daring  both  exposures,  the  dose  was  measured  on 
a  1-ft  rectangular  grid  at  a  height  of  46  in.  above  the  floor.  A  sketch  of  the  salient 
portions  of  the  vent  structure  is  presented  in  Figure  4,  together  with  a  diagram  of 
t^e  two  areas  of  simulated  contamination. 

iJince  the  source  simulation  was  performed  in  rings,  it  is  perhaps  best  to  de¬ 
scribe  the  calculated  infinite-field  dose  in  terms  of  four  annular  areas  of  radii  0  to 
1  ft,  1  to  5  ft,  5  to  25  ft,  and  25  to  *>  ft,  and  to  describe  the  detector  position  in 
terms  of  its  radius  from  the  centerline  of  the  shelter  vent  at  a  height  46  in.  above 
the  shelter  floor.  These  doses  are  further  differentiated  according  to  the  path  by 
which  the  radiation  arrives  at  the  detector. 
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The  following  components  represent  doses  from  the  indicated  radiation  sources 
and  scatterers  to  a  detector  located  at  radius,  r*. 

(0—1)  =  Direct  radiation  from  sources  on  the  vent  cover 

D  (0—1)  =  Radiation  from  sources  on  the  vent  cover  scat¬ 

tered  by  the  walls  of  the  vent 

IX  (0—1)  =  Radiation  from  sources  on  the  vent  cover  scat¬ 

tered  by  the  lower  lip  of  the  vent 

Dj  (0—1)  =  Radiation  from  sources  on  the  vent  cover  scat¬ 

tered  by  the  floor  of  the  shelter 

D0  (1—5)  =  Experimentally-measured  radiation  from  sources 

in  a  contaminated  field  extending  in  a  radius  of  1 
to  5  ft  around  the  vent 

D  (5—25)  =  Experimentally-measured  radiation  from  sources 
in  a  contaminated  field  extending  in  a  radius  of  5 
to  25  ft  around  the  vent 

D  (25— «°)  =  Skyshine  and  vent-scattered  radiation  from  sources 
in  a  contaminated  field  lying  beyond  a  25-ft  radius 
around  the  vent. 

Then  the  total  dose  rate,  (0— «),  will  be  given  by  the  sum  of  these  seven  components. 

Computation  of  the  various  dose  components  requires  data  for  the  direct  and 
scattered  dose  from  each  annulus  and  the  attenuation  afforded  by  different  geometries. 
For  purposes  of  clarity,  the  data  required  are  summarized  at  this  point.  Table  3 
presents  the  direct  and  scattered  dose  components  for  cobalt-60  radiation  from  the 
4  annular  fields  at  different  altitudes  above  the  field.  The  buildup  factors  used  in 
the  computation  were  obtained  from  ref.  5.  The  attenuations  afforded-by  concrete 
(or  steel)  to  cobalt-60  radiation  for  different  geometrical  situations  are  obtained 
from  ref.  3  and  are  shown  in  Figure  3,  p.  5.  We  shall  proceed  with  the  computation 
of  the  dose  components. 

THE  DIRECT  DOSE  FROM  CONTAMINATION  ON  THE  VENT  COVER,  Dd(0—  1) 

It  is  difficult  to  compute  exactly  the  uncollided  dose  reaching  a  detector  located 
at  radius  r  in  the  shelter  from  contamination  existing  on  the  vent  cover,  because 
the  mass  thickness  between  source  and  detector  is  a  rapidly  varying  function  of  the 
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TABLE  3 


DIRECT  RADIATION  INTENSITY  (It/hr)  ABOVE  THE  CENTER  OF  A 
CONTAMINATED  AREA  OF  COBALT -GO,  1  CURIE/FT2 


Contaminated 


Ft  above  Plane 


Amrnlua 

0 

1 

3* 

6 

11.67 

0-1  ft 

— 

31 

4.7  +  0.0 

1.2 

0.3 

1-5  ft 

145 

115 

54+9.8 

20 

6.4 

5-25  ft 

141 

140 

128  +  19 

106 

65 

25-- ft 

212 

211 

211  +  63 

209 

203 

* 

Air-ac altered  intensities  3  ft  above  the  plane  containing  the  annuli  are  given  as 
the  second  figure  in  this  column  (data  based  on  R.  E.  Rexroad  and  M.  A.  Schmoke  ®). 
Infinite -fieW  scattered  radiation  data  are  from  L.  V.  Spencer.2 


radios  of  the  detector  from  the  centerline  of  the  vent.  For  example,  a  detector  cen¬ 
tered  below  the  vent  cover  has  essentially  no  mass  attenuation,  while  a  detector 
several  feet  away  has  an  earth  cover  of  3  ft  blocking  a  portion  of  the  vent  cover 
while  it  stm  maintains  direct  view  of  the  remainder  of  the  vent  cover.  Since  the 
earth -cover  thickness  is  large,  we  may  approximate  the  uncollided  dose  to  a  detec¬ 
tor  by  assamlng  that  all  photons  penetrating  the  earth  are  either  scattered  or  absorbed. 
Thus,  the  dose  rate  to  a  detector  that  is  located  away  from  the  centerline  of  the  vent 
is  a  fund  inn  of  only  that  area  of  the  vent  cover  that  may  be  "seen "  by  the  detector 
<the  source  is  assumed  to  be  uniformly  smeared  on  the  vent  cover)  and  the  mean 
distance  tern  the  detector  io  the  center  of  the  visible  area  of  the  vont  cover.  Then, 
the  ■poolHdsd  dose  reaching  a  detector  from  contamination  existing  on  the  vent 
ocnrer  maj  he  calculated  approximately  by  the  relation: 


Dd  (0-1)  m 


va2qg 


“IlL 

w(0) 


B(x) 


d‘ 


i 


! 

t 


where 

r  =  detector  radius  from  the  vent  axis,  ft 

a  =  radius  of  source  =•  1  ft 

q  =  dose  conversion  factor  =  14.3  R/hr-curie  at  1  ft 

2 

<r  «  source  density  =  1  curie/ ft 

w(r)=  solid  angle  fraction  of  the  portion  of  vent  cover  seen 
by  the  detector  at  radius,  r 

d  =  vertical  distance  from  the  center  of  the  source  area 
to  the  detector  =  14  ft 

B(x)  =  barrier  attenuation  introduce  i  by  0.  25-in,  thick 
vent  cover  =  0. 36. 

Table  4  presents  the  computations  of  the  direct  dose  from  contamination  on 
the  vent  cover. 

TABLE  4 


THE  DOSE  RATE  DUE  TO  CONTAMINATION  ON  COVER 


Radius,  r(ft) 

0 

1 

2 

3 

4 

oj(r) 

.0023 

.0023 

■a 

0 

Dd(0-1)  (R/hr) 

.081 

.080 

H 

.018 

0 

THE  DOSE  FROM  SOURCES  ON  THE  VENT  COVER  SCATTERED  BY  THE  WALLS 
OF  THE  VENT,  Dy  (0-1) 

The  radiation  reflected  from  the  vent  walls  to  a  detector  located  at  a  distance  r 
from  the  centerline  of  the  vent  structure  may  be  computed  by  summing  the  components 
of  the  dose  striking  each  differential  element  of  wall  area  from  each  differential 
element  of  contamination  area  times  the  appropriate  albedos,  over  all  vent  and  wall 
areas.  Since  this  task  would  require  a  computer,  several  simplifying  assumptions 
were  made.  First,  the  source  was  assumed  to  be  a  point  isotropic  source  of  the 
same  strength  as  the  total  cover  contamination  located  at  the  center  of  the  cover. 
Second,  the  cylindrical  vent  was  assumed  to  be  an  octahedral  vent,  and  the  incident 
angle  of  radiation  was  taken  as  the  angle  between  a  perpendicular  line  to  the  center 
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of  each  surface  area  seen  by  the  detector  ami  a  line  from  tin1  center  of  each  area 
through  the  source  locations.  Third,  because  data  were  lacking  for  the  differential 


?  albedo  for  steeL,  smoothed  plots  of  the  differential  albedo  for  concrete  at  1  Mev 

g 

i  were  used.  Figure  5  presents  the  smoothed  curves  of  this  albedo  for  incident 
,  angles  of  74. 5  and  84. 5°.  The  total  dose  received  at  a  detector  from  radiation 
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scattered  from  the  vent  wall  may  tlus  be  written  as  follows  (see  Figure  6  for  the 
angles  and  geometry  involved): 


^  S  q  A.  o(6  ,  ©,  $)  cos  0  B(x)  cos  0 

Dv^l)  -I  - l“V? - 

Q  1 1 


where 

S 

q 

A. 

i 

b 

c 

0 

o 

0 

$ 

ai%.  e.  *) 

B(x) 


source  strength  =  x  curies 

dose  conversion  factor  =  14. 3  R/hr-curie  at  1  ft 

area  of  the  j^1  wall  as  seen  by  tlie  detector 

distance  from  center  of  A.  to  the  source  (ft) 

distance  from  center  of  A.  to  detector  (ft) 

polar  angle  of  incidence  of  gamma  ray  at  center 
of  A.  area 

polar  angle  of  exit  of  gamma  ray  at  center  of  A, 
area  to  detector  1 

azimuthal  angle  of  exit  of  gamma  ray  at  center  of 
jth  area  to  detector 

flux  albedo  -  fraction  of  the  dose  flux  impinging  at 
polar  angle,  0Q,  on  area,  Aj,  that  is  reflected  at 
azimuthal  angle,  *,  and  polar  angle  0,  per  steradian. 
barrier  factor  for  vent  cover  =  0. 36. 


Since  the  approximation  that  all  radiation  from  the  reflecting  surface  strikes 
and  exits  at  the  area's  center  (and  hence  at  this  same  polar  angle)  is  rather  poor  for 
dosimeters  that  can  see  the  entire  vent  wall  (radius  from  0  to  1  ft),  further  compu¬ 
tations  were  performed  by  dividing  the  area  seen  by  the  detector  into  2  and  4  vertical 
sections  and  treating  each  section  as  before.  Table  5  presents  the  dose  rate  contribu¬ 
tions  resulting  from  radiation  reflected  from  the  vent  walls. 


THE  DOSE  FROM  SOURCES  ON  THE  VENT  COVER  SCATTERED  BY  THE  LOWER 
UP  OF  THE  VENT,  D£  (0-1) 

The  dose  was  calculated  for  radiation  from  sources  on  the  vent  cover  striking 

the  lower  lip  portion  of  the  vent  and  scattering  to  the  detector  located  at  radius  r 

7 

from  the  centerline  of  the  vent.  The  procedures  of  LeDoux  and  Chilton  were  used, 
which  are  predicated  upon  a  rectangular-shaped  duct  and  a  point  source.  For 
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TABLE  5 


THE  WALL-SCATTERED  DOSE  RATE  (R/hr) 


Radius,  r(ft) 

0 

1 

, 

2 

4 

i 

6 

8 

10 

Dy  (0—1) 

.0084 

.0069* 

. 005 2 t 

.016 

.012* 

.oiot 

.008 

.0021 

.  0006 

. 00023 

. 00012 

Computed  with  2  vertical  area  sections. 
^Computed  with  4  vertical  area  sections. 


simplicity,  we  assumed  that  the  cylindrical  shelter  vent  was  replaced  by  a  rectangular 
vent  of  similar  size  and  the  area -distributed  source  on  the  vent  cover  was  replaced 
by  a  point  source  of  similar  strength.  Using  LeDoux's  methods,  we  computed  the 
lip-scattered  doses  for  various  detector  positions,  which  are  presented  in  Table  6. 


TABLE  6 

THE  LIP -SCATTERED  DOSE  RATE  (R/hr) 


Radius,  r(ft) 

1 

2 

•  4 

6 

8 

10 

0 

.034 

028 

.015 

.0070 

.0039 

THE  DOSE  FROM  RADIATION  SCATTERED  FROM  THE  FLOOR  OF  THE  SHELTER, 
D{  (0-1) 

The  dose  contribution  from  radiation  scattered  off  the  floor  arriving  at  a  detector 
located  at  radius  r  from  the  centerline  of  the  vent  may  be  approximated  by  summing 
the  direct-dose  rates  at  the  floor  from  each  differential  source  area  on  ihe  vent  cover, 
times  the  appropriate  albedo,  times  each  differential  floor  area,  divided  by  the  dis¬ 
tance  to  the  detector  squared.  As  this  procedure  would  require  a  double  integration 
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over  floor  and  vent  cover  area,  the  following  approximation  was  made.  The  source 
was  assumed  to  lx--  a  p°'nt  source  of  total  strength  equal  to  the  cover  contamination, 
and  the  albedo  was  taken  as  the  albedo  to  the  center  of  the  direct -dose  illuminated 
floor  area.  The  dose  rate  from  floor -scattered  radiation  may  thus  be  approximated 

as: 


Df  (0-1) 


S  q  c*(©0,  ©,  4‘)  cos  0 

F? 


B(*). 


where 

S  =  60 urce  strength  =  it  cur  es 

q  =  dose  conversion  factor  -  14  3  R/hr-curie  at  1  ft 

g  =  distance  in  ft  from  floor  to  vent  cover 

Aj  =  area  of  the  floor  illuminated  by  direct  radiation 

l  =  slant  distance  from  detector  to  center  of  illuminated 

floor  area 

B(x)  -  attenuation  of  vent  cover  -  0. 36 

differential  flux  albedo  per  steradian.  .  . 

The  floor -scattered  dose  was  computed  for  various  detector  locations  and  is 

presented  in  Table  7. 


«(®0.  ©.  *> 


TABLE  7 

THE  FLOOR-SCATTERED  RADIATION  FROM  SOURCES  LOCATED 
ON  THE  VENT  COVER  (R/hr) 


Df  (0-1) 

|  Radius,  r 

(ft) 

0 

1 

2 

4 

6 

8 

'10 

.010 

.0012 

mm 

m 

1 

.011 

.00035 

.014 

.00025 

.016 

. 00015 

.017 

. 00007 

THJE  DOSE  FROM  TWO  ANNULAR  CONTAMINATED  AREAS,  DJI— 5)  AND  D  (5-25) 

©  © 

Figures  7  and  8  show  Shumway's*  experimental  data  from  a  set  of  dosimeters 
arranged  46  in.  above  the  floor  of  the  shelter  and  at  various  distances  from  the  vent 
axle.  Figure  7  represents  measurements,  normalized  to  a  contamination  density  of 
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2 

1  curie  of  cobalt-CO  per  ft  .  obtained  from  an  annular  source  of  1  -ft  inner  radius 
and  5 -ft  outer  radius.  Figure  8  represents  normalized  dose  rates  for  a  source  of 
5-ft  inner  radius  and  25-ft  outer  radius.  In  both  cases,  the  data  showed  no  signifi¬ 
cant  variation  with  azimuth  angle  around  the  vent  axis,  indicating  that  most  of  the 
radiation  reaching  the  detectors  entered  the  shelter  through  the  symmetrical  vent. 
This  radiation  can  be  scattered  either  by  the  air  above  the  vent  or  by  the  vent  struc¬ 
ture  itself.  We  have  calculated  these  two  types  of  contributions  to  obtain  their 
relative  magnitudes  and  to  assist  in  making  the  necessary  extrapolation  to  infinite- 
contaminated-field  conditions.  We  shall  discuss  these  calculations  for  skyshine 
and  vent -scattered  radiation  in  the  following  paragraphs. 
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Figure  7.  Observed  Dose  Rates  Compared  with  Calculated  Scatter  from  Con¬ 
taminated  Annuli  1-  to  5-Ft  Inner  and  Outer  Radii,  4G  Inches  Above  Floor 
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Figure  8.  Observed  Done  Kates  Compared  with  Calculated  Scatter  from  Con¬ 
taminated  Annuli  5-  to  25-Ft  Inner  and  Outer  Radii,  46  Inches  Above  Floor 

Skyahine 

Several  sources  of  information  may  be  used  to  estimate  the  skyshine  contribu- 

Q 

tion.  Schumchyk  and  Tiller  have  reported  that  measurements  of  scattered  doses 
from  point  sources  of  cobalt -60  lying  on  the  ground  at  various  distances  from  a 

O 

foxhole  agree  reasonably  well  with  theoretical  predictions  by  Spencer.  In  particular, 
they  showed  that  for  sources  between  2  and  1G  ft  from  the  detector,  the  quantity 
Q(w)  =  2*r  Ij.(w)  dependt;  only  on  the  solid  angle  fraction  u;  seen  by  the  detector. 

Here  r  is  the  source  distance  from  the  axis  of  the  foxhole,  and  Ij.  is  the  observed 
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dose  rate  in  R/hr  per  curie.  They  also  found  that  for  small  solid  angles  («  <  0. 2) 
the  variation  of  intensity  with  solid  angle  oould  be  expressed  as 


where  Soglo  is  the  scattered  radiation  present  at  the  top  of  the  hole,  calculated 
from  theoretical  Infinite-medium  buildup  data.  This  same  ratio  was  reported  by 
Clarke,  et  &1. ,  to  represent  Monte  Carlo  calculations  for  predicting  the  scattered 
radiation  entering  an  open  basement  surrounded  by  a  cobalt-60  contamination  field. 

Now,  it  can  be  readily  shown  that  the  integrated  skyahine  dose  rate  on  the  axis 
of  a  hole  surrounded  by  a  contaminated  annulus  with  radii  a  and  b  is 

b  b 

Ig(w)  =  y  (2irrodr)  ^(w)  ■  »  J  Q(u)  dr 
a  a 

"  °<b  -  a)  Q(w) 

since  for  small  radii,  Q  was  shown  to  be  independent  of  r.  Here  <r  roproonnts  the 
source  density  on  the  ground.  Inserting  data  from  ref.  6,  the  intensity  on  the 
is 

I  (oj)  =  57  wo(b  -  a)  mR/hr, 
for  u  <  o.  2,  small  radii,  and  a  in  curles/ft^. 

For  the  two  annuli  used  here,  the  predicted  skyahine  46  in.  off  the  floor  attenuated 
by  a  factor  of  0. 36  by  the  steel  vent  and  lid  is 


1-5  ft  annulus:  I 

s 

5-25  ft  annulus:  I 


0. 31  mR/hr 
1. 55  mR/hr. 


Vent  Scatter 


The  radiation  scattered  into  the  shelter  by  the  vent  was  estimated  with  the  aid 

a 

of  Monte  Carlo  calculations  by  Raso  giving  the  reflection  and  transmission  of  1-Mev 
gamma  radiation  striking  ooncrete.  These  data  present  foe  differential  albedoe  and 
transmissions  as  functions  of  both  polar  and  azimuth  exit  angles  for  several  incident 
polar  angles. 
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To  simplify  the  problem,  it  was  assumed  that:  (1)  all  incident  radiation  arriving 
at  the  vent  travelled  parallel  to  the  ground,  (2)  for  radiation  to  reach  a  detector  in 
the  shelter  and  on  the  cylinder  axis,  the  exit  polar  angle  averaged  85°  and  the  exit 
azimuth  was  always  90°,  (3)  the  scattered  radiation  was  composed  partly  of  gamma 
rays  passing  through  the  cylinder  and  emerging  in  the  downward  direction  and  partly 
of  radiation  passing  through  one  side  of  the  cylinder  and  scattering  backward  and 
downward  from  the  opposite  wall. 

The  equation  used  was  one  presented  in  the  Addendum  to  ref.  C,  giving  the  per 
unit  incident  dose  rate  at  a  point  of  specified  distance  and  direction  from  a  scatterer. 
The  incident  radiation  from  the  contaminated  annuli  was  considered  to  be  almost 
entirely  direct;  its  intensity  is 

D  ■  2irqo  E^lAj h2  +  a2)  -  Ej  (/I\Jb2+  b2), 

where  h  is  the  detector  height,  a  and  b  are  inner  and  outer  annulus  radii,  is  the 

-3 

air.  absorption  coefficient  =  2. 24  x  10  per  ft,  and  the  remaining  terms  are  those 
used  earlier.  Results  of  the  application  of  this  equation  have  been  given  in  Table  3, 

p.  11. 

Dose  rates  due  to  vent-scattered  radiation  at  detectors  directly  beneath  the  vent 
and  46  in.  above  the  floor  were  calculated  to  be  17  mR/hr  and  21  mR/hr  for  the  1  to 
5-ft  annulus  and  the  5  to  25-ft  annulus  respectively.  These  rates  may  be  compared 
with  the  observed  value  of  16  mR/hr  found  for  both  annuli.  In  addition,  for  radiation 
from  the  larger  annulus  scattering  to  a  detector  12  ft  above  the  shelter  floor,  the 
observed  30  mR/hr  may  be  compared  with  the  computed  46  mR/hr.  Differences 
between  measurements  and  calculations  are  in  part  attributable  to  the  use  of  albedos 
for  1  Mev  on  concrete  instead  of  the  actual  case  of  cobalt-60  (1. 25  Mev)  on  iron  and 
in  part  attributable  to  the  assumption  of  a  semi-infinite  backscattering  medium 
rather  than  the  actual  l/3-mfp-thick  vent  wall. 

Radial  Distribution  of  Scattered  Dose 

The  previous  calculations  were  used  to  estimate  the  dose  rate  scattered  to  a 
point  on  the  axle  of  the  cylindrical  vent.  As  an  approximation,  it  was  assumed  that 
the  scattered  radiation  from  either  the  air  or  the  vent  was  distributed  isotropically 


and,  hence,  that  the  dose  rate  would  be  proportional  to  the  fraction  of  the  vent  area 
that  could  be  seen  by  the  detector.  The  problem  is  simplified  to  one  of  determining 
the  fractional  area  common  to  two  overlapping  circles.  These  circles  can  be 
visualized  by  imagining  the  appearance  of  the  vent  when  looked  at  from  below  and  to 
one  side,  as  illustrated  in  Figure  9.  Then,  the  fractional  area  is: 

f  |^cos_1y -y"VT-y2| 

where  y  =  ~  =  j  -  1),  permitting  the  intensity  (as  a  fraction  of  axis  value)  to 
be  calculated  as  a  function  of  distance,  r,  from  the  axis. 


Figure  9.  Geometry  of  the  Shelter  Vent 
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The  results  of  these  calculations  are  shown  in  Figures  7  and  8,  pp.  1 8  and  19. 
There  seems  to  be  satisfactory  agreement  between  the  observed  and  calculated 
intensities,  except  at  points  at  a  considerable  distance  from  the  vent  axis.  Undoubt¬ 
edly,  the  increased  dose  rates  observed  here  are  produced  by  further  scattering 

from  the  vent  walls. 

THE  SKYSHINE  AND  VENT -SCATTERED  DOSE  FROM  THE  CONTAMINATED 
AREA  BEYOND  25  FT,  De  (25—  ») 

The  contribution  to  the  shelter  dose  from  radiation  sources  beyond  the  experi¬ 
mentally  tested  area  must  be  calculated  and  added  to  the  observed  intensities.  The 
dose  rate  scattered  into  the  shelter  by  the  vent  is  ca3y  to  compute;  the  scattered 
dose  is  proportional  to  the  direct  radiation  incident  on  the  vent.  A  vent-scattered 
dose  rate  of  32  mR/hr  from  the  far -field  beyond  25  ft  is  thus  obtained  from  the 
scattered-dose  rates  due  to  the  annuli. 

Skyshine  is  obtained  in  a  different  way.  Table  3  (p.  11)  gives  the  total  skyshine 
for  a  contaminated  field  beyond  25  ft  as  63  R/hr;  this  was  obtained  by  subtracting 
the  annular  rates  from  the  theoretical  scatter  due  to  an  infinite  field  (ref.  3,  Figs. 
28.3  and  B-15).  In  penetrating  to  the  detector  in  the  shelter,  this  skyshine  is 
attenuated  (0. 36)  by  the  thickness  of  the  steel  vent  and  by  the  solid-angle  fraction 
(0. 004)  subtended  by  the  vent  opening.  This  solid-angle  fraction  produces  a 
geometrical  attenuation  of  0. 001  shown  as  an  extrapolation  of  Figure  B.  37  of  ref.  3. 
The  net  result  is  a  calculated  skyshine  of  23  mR/hr  due  to  contamination  beyond 
26  ft. 

The  combined  dose  rate  of  55  mR/hr  represents  the  expected  intensity  on  the 
vent  axis  46  in.  above  the  floor.  As  an  approximation,  it  is  expected  that  the  dis¬ 
tribution  of  intensity  with  radius  from  the  axis  should  follow  the  same  pattern  as 

that  for  the  annuli. 

TOTAL  DOSE  RATE,  Dt  (0—«) 

The  total  dose  rate  in  the  shelter  due  to  the  radiation  from  an  infinite  field  of 
2 

cobalt -60  (1  curie/ft  )  penetrating  the  shelter  vent  can  now  be  estimated  by  sum¬ 
ming  the  various  contributions.  These  are  summarized  in  Table  8. 


TABLE  8 


DOSE  RATES  (mR/hr)  46  IN.  ABOVE  THE  SHELTER  FLOOR  DUE  TO 
1  CURIE/FT2  COBALT -60  ON  SURFACE* 


Radius,  r(ft) 


0 

1 

2 

4 

6 

8 

10 

Dd  (0-1) 

81 

80 

56 

0 

0 

0 

0 

»¥  (o-i) 

5 

10 

8 

2 

0.6 

0.2 

0.1 

Dj  ((►-!) 

0 

0 

34 

28 

15 

7 

4 

Df  (0-1) 

1.2 

0.8 

0.7 

0.35 

0. 25 

0.15 

0. 07 

De  (1-5)  t 

15.4 

14.5 

13 

9.5 

5.8 

2.8 

1.0 

De  (5-25)  t 

17 

13 

10 

5.8 

3.4 

2.0 

1.1 

Dg  (25— «) 

55 

55 

43 

20 

2.5 

0 

0 

Dt  (0— ») 

175 

173 

165 

66 

27 

12 

6 

*Surface  dose  rate  is  500  R/hr  at  3  ft, 

^De  values  are  means  of  experimental  data. 


SHELTER  ENTRANCE  WAY 

The  NRDL  underground  shelter  entranceway  is  an  L-shaped  passage  designed 
to  have  no  line-of-sight  leakage  path  for  radiation.  The  entranceway  is  composed 
of  a  7-ft-diameter  tube  inclined  at  a  45°  angle  to  the  earth's  surface,  penetrating 
to  a  depth  of  approximately  12  ft  where  it  is  connected  to  a  7-ft-diameter  tube  run¬ 
ning  horizontally  for  approximately  12  ft  to  the  shelter.  The  outer  portal  of  the 
entranceway  is  covered  with  a  3/16 -in.  -thick  cover  plate.  Figure  10  illustrates  the 
pertinent  features  of  the  entranceway. 

Shumway*  of  NRDL  exposed  the  entranccway  to  a  rectangular  contaminated  field 
of  cobalt -60.  Two  experiments  were  performed  with  this  source  array.  In  the  first 
experiment,  contamination  was  simulated  over  the  entire  area  including  the  entrance¬ 
way  cover.  The  second  experiment  was  a  duplicate  of  the  first,  except  that  the 
entranceway  cover  was  left  uncontaminated.  This  pair  of  experiments  provides  an 
opportunity  to  test  the  methods  for  computing  radiation  penetration  of  right-angle 
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7*  DIAMETER 


Figure  10.  The  Shelter  Entranceway 
7  12 

ducts  suggested  by  LeDoux,  et  al. ,  and  Eisenhauer.  Shumway's  observations 
along  the  centerline  of  the  horizontal  tube  are  presented  in  Table  9  as  a  function  of 
distance  from  the  junction  of  the  centerlines  of  the  two  legs  of  the  duct.  The  differ¬ 
ences  should  be  attributable  to  contamination  on  the  entranceway  cover. 

12 

Eisenhauer,  in  his  technical  note  on  the  scattering  of  radiation  in  ducts,  pre¬ 
sents  experimentally -measured  data  on  3  ducts  of  rectangular  cross  section  contain 
inf  a  right- angle  bend.  He  shows  from  single-scatter  theory  that  the  dose  in  the 
seoond  leg  should  be  expressible  by  the  equation 


where 


D  =  dose  rate  along  the  centerline  of  the  horizontal  leg 

D  =  dose  rate  at  the  junction  of  the  centerlines  of  the  first 
°  and  second  leg 

i o  =  the  sol  id -angle  fraction  subtended  at  the  detector  by 
the  junction. 

When  the  solid  angle  was  defined  in  terms  of  the  duct  cross  section  and  the  detector 
distance  from  the  near  side  of  the  intersection,  his  experimentally -obtained  data 
were  found  to  agree  with  this  equation  for  distances  of  several  diameters  from  the 
junction. 


! 


f 

1 


TABLE  9 

OBSERVED  DOSE  RATES  IN  THE  HORIZONTAL  ENTRANCE  WAY  DUE  TO 
AREAS  CONTAMINATED  WITH  1  CURIE/FT2  COBALT -60* 


Distance  from 
Centerline  Junction 
(ft) 

Entry  Cover 
Contaminated 
(R/hr) 

Entry  Cover 
Clear  (R/hr) 

Difference 

(R/hr) 

0 

2.60 

0.  480 

2.12 

2.8 

2.  70 

0. 510 

2. 19 

3.5 

1.90 

0.  460 

1.44 

5.5 

0. 330 

0.114 

0.  216 

7.5 

0. 044 

0.025 

9 

0. 022 

0.  0073 

1' 

12 

0. 0091 

0.  0023 

1 

15 

0. 0053 

0. 0019 

0. 0034 

Surface  dose  rate  is  500  R/hr  at  3  ft. 


7 

A  second  approach  to  the  problem  is  presented  by  LeDoux,  who  divides  the 
junction  of  the  two  ducts  into  several  suitable  areas  and  computes  the  dose  contribu¬ 
tion  of  each  area  by  a  straightforward  albedo  approach.  Correction  terms  aro 
applied  to  account  for  the  corner  "lip  effect.  "  Because  differential  albedo  data 
were  not  available  at  the  time,  LeDoux  used  the  results  of  Monte  Carlo  calculations 
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by  Berger  and  Raso,*^  who  reported  albedos  of  surfaces  irradiated  at  various  incident 
polar  angles.  LeDoux  assumed  that  the  backscattcred  currents  were  distributed 
isotropically  and  calculated  differential  albedos  accordingly.  The  agreement  vith 
Eisenhauer's  experimental  flata  was  good.  Raso,  however,  has  since  calculated 
differential  dose  albedos  of  1-Mev  photons  incident  on  concrete,  and  we  therefore 
attempted  to  modify  LeDoux's  work  with  this  more  recent  information.  A  plot  of 
Raso ’a  data  on  differential  dose  albedos  applicable  to  this  calculation  is  shown  in 
Figure  11,  along  with  curves  representing  LeDoux's  assumed  values.  The  particular 
directions  were  chosen  because  they  cover  the  four  scattering  surfaces  at  the  elbow 
of  a  square  duct. 

For  computational  purposes,  the  cylindrical  entranceway  of  the  NRDL  shelter 
was  replaced  by  a  square  entranceway  of  equal  cross-sectional  area,  since  both 
methods  described  above  were  based  on  rectangular  ducts.  The  results  of  these 
calculations  are  shown  in  Table  10  in  terms  of  the  ratio  of  observed-to-calculated 
reductions  in  intensity  relative  to  the  radiation  intensity  incident  on  the  junction 
midpoint.  (The  measured  intensity  of  the  radiation  incident  on  the  junction  agreed 
closely  with  the  valuo  calculated  for  the  contaminated  entry  cover. ) 


TABLE  10 

COMPARISON  BETWEEN  CALCULATED  AND  OBSERVED  INTENSITIES 
ALONG  THE  HORIZONTAL  ENTRANCE  WAY 


All  three  computations  produced  larger  scattered-dose  rates  than  were  actually 
observed.  The  agreement  between  Eisenhauer's  and  LeDoux's  isotropic-computed 
doee  rates  was  to  be  expected,  since  both  had  been  checked  against  Eisenhauer's 


DIFFERENTIAL  ALBEDO  a(0o,e,  <$ )  PER  STERADIAN 


Figure  11.  Differential  Dose  Albedos  of  1  Mev  on  Concrete,  from 
Several  Particular  Directions 


experimental  data.  The  difference  between  the  two  IxjDoux  calculations  is  due  to 
differences  in  the  assumed  albedos,  which  in  part  may  be  ascribed  to  the  use  of 
1  Mev  rather  than  1.  25  Mev  Monte  Carlo  data.  But  there  is  clearly  a  factor  of 
from  2  to  3  between  the  NRDL  shelter  measurements  and  those  originally  made  on 
the  concrete  ducts. 
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The  extra  shielding  demonstrated  by  the  NRDL  design  may  re  it  from  the 
corrugated  shape  of  the  tube  comprising  the  entranceway.  The  effect  of  these  cor¬ 
rugations  would  be  to  make  the  incident  angle  of  the  primary  radiation  more  nearly 
normal,  thereby  reducing  the  albedo  from  the  higher  values  characteristic  of  grazing 
radiation.  The  fact  that  part  of  the  surface  is  shielded  by  these  corrugations  does 
not  affect  the  albedo,  since  the  total  flux  striking  the  surface  is  not  changed  by  sur¬ 
face  roughness. 

It  is  also  of  interest  to  estimate  the  effect  of  far-fic.ld  skyshine  radiation  pene¬ 
trating  the  shelter  entranceway.  An  approximate  method  of  estimating  the  effects 
of  far-field  skyshine  radiation  has  been  previously  discussed  (p.  23).  Skyshine 
radiation  down  a  hole  was  shown  to  follow  the  relationship, 

Vj  =  -^p-^  for  a)  <  0.  3, 

where 

I  =  the  skyshine-dose  rate 

I  =  the  scattered-dose  rate  at  the  entrance  to  the  hole 

o 

w  =  solid-angle  fraction  of  the  mouth  of  the  hole  from  the  detector 

B(x)  =  barrier  factor  introduced  by  the  hole  cover  =  0.  43  for  a  3/16- 
in.  -thick  iron  cover. 

The  quantity,  IQ,  may  be  approximated  by  assuming  that  the  25  x  50-ft  rectangular 
source  area  is  replaced  by  an  annular  area  of  the  same  size,  centered  on  the  entrance¬ 
way,  and  calculating  the  scattered-dose  rate  to  be  expected  from  sources  lying 

5 

beyond  that  area.  Using  the  buildup  factors  measured  by  Rexroad,  with  the  source 
at  ground  level  and  the  detector  at  1  ft  height  as  before,  IQ  is  found  to  be  equal  to 
68. 3  R/hr. 

Table  11  presents  the  calculated  values  of  far-field  skyshine  together  with  the 
experimentally-measured  values  of  dose  contribution  from  a  25  x  50-ft  rectangular 
source  area  with  the  entranceway  cover  cleared,  the  dose  contribution  from  a  con¬ 
taminated  entranceway  cover,  and  the  total  of  all  contributions. 
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TABLE  11 


TOTAL  DOSE  RATE  (Il/hr)  IN  THE  SHELTER  ENTRANCEWAY  FOR  A 
UNIFORMLY-CONTAMINATED  INFINITE  FIELD  OF 
1  CURIE/FT2  COBALT -60* 


Centerline 

Distance  From  Entranceway  Cover,  Ft 

7 

10 

13 

16.25 

Far-Field  Skyshine 

1.1 

0.6 

0.37 

0.  24 

Inner  Rectangular  Area 

4.9 

2.0 

0.95 

0.  48 

Entranceway  Cover 

13.1 

5.2 

3.25 

2.12 

Total  Dose 

19.1 

7.8 

4.  57 

2.  84 

$ 

Surface  dose  rate  is  500  R/hr  at  3  ft. 


GROUND  PENETRATION 

The  NRDL  underground  shelter  resembles  a  quonset  hut  approximately  49  ft 
long,  with  an  arch  of  approximately  12.  5-ft  radius.  A  clay-like  substance  with  a 
density  of  1. 68  when  compacted,  1. 37  when  uncompacted,  and  1.  46  when  wet  was 
used  to  cover  the  shelter.  The  minimum  earth  cover  at  the  center  of  the  arch 
varied  from  3  ft  8  in.  near  the  enti  anceway  of  the  quonset  hut  to  2  ft  10  in.  near  the 
rear  of  the  structure.  The  earth  was  allowed  to  settle  naturally  and  desiccate.  For 
computational  purposes,  a  density  of  1. 46  will  be  used. 

Experimental  measurements  of  the  attenuation  afforded  by  this  earth  cover 
against  cobalt -60  radiation  were  performed  and  reported  by  Shumway.*  The  experi¬ 
ment  consisted  of  centering  a  simulated,  uniformly-contaminated  30  x  30-ft  field 
over  the  shelter  with  one  edge  in  the  plane  of  the  shelter's  front  wall  (entranceway 
wall),  as  shown  in  Figure  12. 

Because  a  portion  of  the  array  could  shine  into  the  entranceway,  only  data  due 
to  sources  kept  removed  from  the  entranceway  were  used  for  comparison  of  com¬ 
putational  and  experimental  results.  Dose  rates  of  1  to  2  mR/hr  with  an  uncertainty 
of  ±0.4  mR/hr  were  measured  from  a  field  of  1  curie/ft  source  density.  The  depth 
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of  earth  cover  varied  over  the  experimental  area  from  a  maximum  of  3  ft  8  in.  to  a 
minimum  of  3  ft  5  in.  over  the  center  of  the  arch.  Thus,  for  computational  purposes, 
an  average  value  of  3  ft  6  in.  will  be  used.  This  value,  together  with  the  best  estimate 
of  earth  density  (p  =  1. 46),  provides  an  estimated  average  value  of  earth  cover  of 
320  lbs/ft2  at  the  arch  center. 

Three  methods  of  computing  the  attenuation  of  an  arch-type  shelter  are  avail  able 

13 

in  the  published  literature.  One  method  is  by  LeDoux  of  the  U.  S.  Naval  Civil 
Engineering  Laboratory.  LeDoux  proposes  that  the  attenuation  afforded  by  an  under¬ 
ground  shelter  is  equivalent  to  the  attenuation  of  an  earthen  slab  of  thickness  equal 
to  the  shelter  minimum  roof  thickness  times  a  geometry  factor.  The  geometry 
factor  is  based  on  an  analysis  of  radiation  penetrating  the  shelter  from  discrete 
sources  lying  on  the  ground  above,  using  an  exponential  approximation  for  the  build¬ 
up  factor.  For  a  horizontal  semicircular  cylindrical  shelter,  his  results  lead  to  a 
geometry  factor  of 

I  a  +  _10 

6  2  +  h/t>' 

where  h  is  the  radius  of  the  shelter  arch,  and  t  is  the  minimum  earth -cover  distance. 
This  factor,  effectively  independent  of  energy,  applies  to  the  intensity  on  the  floor 
of  the  shelter. 

15 

Malich  and  Beach  state  that  the  dose  rates  at  the  centerline  of  an  arched 
shelter  3  ft  above  the  floor  are  well  represented  by  90%  of  the  dose  rate  that  would 
be  received  in  a  shelter  with  a  slab  roof  of  thickness  equal  to  the  minimum  earth 
cover  of  an  arched  shelter.  Malich  and  Beach  attribute  this  to  the  fact  that  over 
90%  of  the  dose  rate  within  a  shelter  comes  from  a  relatively  small  circle  on  the 
ground  directly  above  the  detector.  The  radius  of  this  circle  ranges  from  2. 6  ft  at 
a  depth  of  3  ft  to  4. 8  ft  at  a  depth  of  10  ft.  Thus,  we  may  assume  that  the  experi¬ 
mental  30  x  30-ft  source  area  adequately  represented  an  infinite  field  for  ground 
penetration  measurements. 

3 

A  third  method  by  Spencer  pi’ovides  detailed  information  for  both  cobalt -GO 
and  1. 12-hr  fallout  radiation  on  the  angular  distribution  of  the  radiation  dose  pene¬ 
trating  various  thicknesses  of  earth,  thus  allowing  the  summation  of  radiation  com¬ 
ponents  from  each  individual  increment  of  source  area.  Since  the  NRDL  shelter 
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was  symmetric  about  the  vertical  plane  extending  from  the  floor  through  the.  shelter 

arch,  differential  strips  of  1  -ft  wide  contamination  extending  the  length  of  the  shelter  i 

were  used  in  computing  the  dose  rate  on  the  centerline  3  ft  above  the  shelter  floor.  j 

i 

Spencer  (ref.  3,  Fig.  B-13)  also  provides  information  on  the  attenuation  provided 

by  a  slab  covered  with  cobalt-60:  for  320  psf,  the  reduction  factor  relative  to  the 

“5 

dose  rate  at  3  ft  above  the  contaminated  plane  is  4.  5x10  .  This  will  reduce  an 

external  field  of  500  R/hr  at  3  ft  to  about  22  mlt/hr.  Table  12  compares  the  dose 
rates  for  the  cylindrical  shelter  calculated  using  these  three  methods  with  the  dose 

rates  obtained  experimentally. 


TABLE  12 

GROUND  PENETRATION  OF  COBALT -60  RADIATION  FOR 
500  R/hr  EXTERNAL  FIELD  (1  CURIE/FT2) 


Experimental 

(Shumway) 

1  Calculated 

LeDoux 

Malich  &  Beach 

Spencer 

Geometry  Factor 

— 

0.47 

0. 90 

— 

mR/hr  in  Shelter 

2 

10 

20 

16 

Reduction  Factor 

4  x  10~6 

2  X  10'5 

4  x  10*5 

3  x  10~5 

The  agreement  between  experiment  and  computation  is  poor.  Even  if  we  assume 
that  the  sofl  density  is  the  maximum  compacted  value  producing  a  minimum  over¬ 
head  mass  of  370  psf,  attenuation  is  increased  by  only  a  factor  of  3.  To  try  to 
resolve  the  disagreement,  we  analyzed  data  obtained  by  Strope11  on  a  buried  quonset- 
type  shelter  similar  to  the  NRDL  shelter,  which  was  subjected  to  actual  fallout  from 
a  nuclear  weapon  test.  The  results  of  these  computations,  shown  in  Table  13,  are 
In  fair  agreement  with  the  experimental  data.  Thus,  the  discrepancy  between  com¬ 
puted  and  measured  attenuation  for  the  NRDL  shelter  may  probably  be  attributed  to 
an  error  in  the  measurement  of  earth-cover  thickness  or  density. 

ft  Is  te  be  noted  that  the  reduction  factors  for  ground  penetration  of  fallout 
radiation  are  appreciably  greater  than  those  for  penetration  of  cobalt-60  radiation 
and  that,  aider  attack  conditions,  the  radiation  penetrating  the  interior  of  the 
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efielter  through  the  roof  may  not  be  negligible  compared  to  that  coming  through  the 
vent  and  the  entranceway. 


TABLE  13 

GROUND  PENETRATION  OF  FALLOUT  RADIATION  FOR 
2. 2  R/hr  EXTERNAL  FIELD 


CHAPTER  4 


CONCLUSIONS 

Table  14  summarizes  the  estimates  of  radiation  levels  expected  to  prevail  in  the 

underground  shelter  if  the  ground  above  were  covered  with  an  infinite  plane  of  cobalt- 

2 

60  of  density  1  curie/ft  ,  corresponding  to  500  R/hr  measured  3  ft  above  the  plane. 


TABLE  14 

INTERIOR  DOSE  RATES  DUE  TO  COBALT -60  AND  FALLOUT  FIELDS 
OF  500  R/hr  AT  3  FT  ABOVE  SURFACE 


Mode  of  Penetration  ! 

Point  of  Measurement 

Intensity  (mR/hr) 

(3  ft  above  floor) 

Cobalt -GO 

Fallout 

Through  roof 

Throughout  shelter 

2 

65 

Directly  below  vent 

175 

225 

Through  vent 

6  ft  away 

27 

40 

10  ft  away 

6 

9 

Through  entranceway 

At  lower  door  into  shelter 

22 

33 

(lid  contaminated) 

6  ft  into  shelter 

5 

8 

Through  entranceway 

At  lower  door  into  shelter 

7 

10 

(lid  cleared) 

6  ft  into  shelter 

2 

3 

By  applying  the  appropriate  factors  from  Table  2  (p.  7),  we  can  arrive  at  the 
expected  intensities  to  be  found  in  the  shelter  when  equal  intensity  fallout  radiation 
(500  B/hr  3  ft  above  the  surface)  is  generated  by  close-in  fallout  1  hr  old. 


It  may  be  concluded  that  over  most  of  its  floor  area  the  shelter  offers  protection 
factors  of  100, 000  against  fallout  radiation  and  factors  greater  than  10,  000  every¬ 
where  except  within  6  ft  of  the  vent  centerline.  The  two  weak  areas  are  under  the 
vent  and  near  the  entranceway.  If  contamination  can  be  cleared  from  the  lid  of  the 
entranceway,  the  levels  near  the  door  can  be  halved. 
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